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I review the linear and second-order perturbation theory in dark energy models with explicit in- 
teraction to matter in view of application to A'^-body simulations and non-linear phenomena. Several 
new or generalized results are obtained: the general equations for the linear perturbation growth; 
an analytical expression for the bias induced by a species-dependent interaction; the Yukawa cor- 
rection to the gravitational potential due to dark energy interaction; the second-order perturbation 
equations in coupled dark energy and their Newtonian limit. I also show that a density-dependent 
effective dark energy mass arises if the dark energy coupling is varying. 



INTRODUCTION 



O ■ 

■ 

Dark energy is defined as a fiuid distributed almost homogeneously and whose potential energy dominates over 
^ ' kinetic energy. As such, it can be observed mainly through large scale eflFects as those relating to the cosmic expan- 
O ' sion history and its linear fluctuations. Indeed, the most weighty evidences in favor of dark energy come from the 
• acceleration of the universe as seen in the Hubble diagram of the supernovae la Q and on the angular size of the 

1 \ acoustic horizon on the cosmic microwave background [21 . Indications concerning the growth of linear fluctuations, 
I e.g. via the ISW effect 0, 0| are still very tentative, although the prospects from e.g. weak lensing Q and Lyman-a 

clustering 0, 01 appear promising. 
. However, all these observables depend ultimately on dark energy only through the expansion history H{z) and the 
I> ' matter linear growth function D{z), where z is the cosmological redshift. For instance, luminosity distance in flat 
! space is defined as 



_ where 

. wViPrp n — (^ 4- 

to 



d.iz)- (1 + ^)1^ (1) 



H{z) = Ho [f7„a-3 + 17^a-3[i+^('^)l + (1 - f7„ - n^)a-^] , (2) 



where a — {1 + z) ^ is the scale factor, fli denotes the density at the present time of the i-th species and 

p3 ■ W{a = log a) = - / w4a')da', (3) 

' Jo 



k> , w^{z) being the equation of state of dark energy. It is clear that at any given redshift there will be different w^{z) 
^ ■ that give indistinguishable dL{z) and that the degree of degeneration will increase with redshift. Similar integrals of 
■ H{z) enter the definitions of angular-diameter distance and age, that will therefore be subject to the same ambiguity. 
The linear growth function is a second observable quantity, in general independent of H{z). In a matter dominated 
epoch with present density parameter and on sub-horizon scales it is given by the solution of the perturbation 
equation 

1 / , . ^1 3f2„iO -3, 



D" + ^[^ + ^)D'-^a-'D = (4) 

(valid for uncoupled dark energy) where the prime denotes derivation with respect to a = log a and where we introduce 
the conformal Hubble function H — da/dr, r being the conformal time. The growth D is therefore an independent 
probe of dark energy: two models that give an identical H{z) will in general be distinguished by different D{z). 

In principle the degeneracy can be broken by a large number of observations at different z of H{z) and/or D{z). 
However, this is hardly feasible, since real data is confined to "small" (z < 5) or very large redshifts {z w 1100). 
Moreover, in most models the dark energy component become subdominant at z 1, so that both H{z) and D{z) 
becomes rapidly insensitive to w^{z) (see for a detailed discussion on the practical observability of w^{z) at large 
z). 
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It would be desirable therefore to add new observables, as the evolution of the perturbations in the dark energy 
field itself, the second order growth function or the full non-linear properties as obtained through A'^-body simulations 
(see e.g. |aUiilll|)- In this paper we review the general linear and non-linear perturbation equations in dark energy 
models in order to provide the basic material for the study of these additional observable quantities. Our dark energy 
model is quite general: a scalar field with a generic potential and an explicit varying species-dependent coupling to 
baryons and dark matter. This cover most scalar field models presented in literature, with the notable exception of 
models with non-standard kinetic terms [12]. 

In a subsequent paper we employ the formaHsm derived here to evaluate the large scale skewness in coupled dark 
energy. 



II. COUPLED DARK ENERGY 



Our dark energy model is characterized by a general potential V{(j)) and general couplings Cj (6) to matter. This 
class of models include those motivated by string theory proposed in 0, Q, 0, 0, 0, O, llfll l2(l Ell l23 |. The 
conservation equations with interacting terms for the field , cold dark matter (c), baryons (6 ) are: 

Ttc).;, = -Ccimc)^;. , (5) 

'^(b)^;t. = -Cfc(0)r(fc)0;^ , (6) 

= [Cc(0)r(e)-fa(0)T(,)]0;,, (7) 

where the coupling functions Cb,c{4') depend on the species, as first proposed in Radiation (subscript 7) remains 
uncoupled because it is traceless (models with coupling to the electromagnetic field [3, or neutrinos [i^l have 
also been proposed). The standard Einstein equations are assumed to hold. This coupling form is derived, through a 
conformal transformation, from a Brans-Dicke gravity with species-dependent interaction [jS. l2^. There are strong 
limits on the baryon coupling [2^ and relatively looser ones on the dark matter coupling from either astrophysics 
or cosmology fs^]; here, however, we disregard any constraint and leave the couplings as free functions. In a 
FRW metric the equations Q plus the Friedmann equation read: 

Pc + 3Hpc = -\/2jiK.I3cPc4>, 

Pb + 3Hpb = -y/2/SKPbpb<P, (8) 
Pj + AHp^ = 0, 

3iJ^ = (/O7 + Pc+ Pb+ P^) , 

where = SttG, /3c — Cc\J Pb = ^b\J 5^1 H — ^l^^- The matter conservation equations can be integrated out: 

Pc,fc = p(0)c,ha"^cxp[-y^K j (3c,b{4>)d4>\- (9) 

This shows one of the basic properties of dark energy interactions: matter density appears to be non-conserved. For 
as concerns the potential V{(j)) here we write in all generality 

V{(t)) = Ae~''N/273f/(«<#> , (10) 

where /i is a dimensionless constant. The exponential case studied in [T^. Il3 | corresponds therefore to / = 1 , a 
constant potential as in to Ai = 0, and the power law V - 0"" to /(</>) = ^(3/2)nlog(/)/(K/^0). We will use also 
the definitions 
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where 

h = + (13) 

h-f!-^'^- (14) 

K/i d(p 

The simple exponential case reduces then to / = /i = /2 = 1. We'll need later on higher derivatives of V so we give 
the general rule: 

v^;^ = i-ir (^) ^ ^>"/„^ , 

where fo = l and 

/„ = fn-lh - fn-l/(.\ll^f^) ■ (15) 

The system ||Sl is best studied in the new variables [T^l3ll| 

and the time variable a = log a. Then we obtain 

■ - 1 ) X - Ai/iy' + /?c(l -x^-v'^-v'' - z^) + I3bv^, 



y' = fJ-fixy + y [2 



z 

z' = -^(l-3.x2 + 3y2-z2), (17) 
v' = {S(3bx - 3x'^ + 3y^ - z'^) . 

The CDM energy density parameter is obviously flc ^ I — ~ y"^ ~ — while we also have fl^ — x^ + y"^ , fly — z^ 
and rib = v^- The system is subject to the condition x^ + + v"^ + z"^ < 1. To close the system one needs also the 
relation fi{y,H) and the Friedman equation 

^ = -\{3 + 3x'-3y' + z'). (18) 

In the perturbation calculations we will use the conformal Hubble function, H — aH, and we will need the following 
relation 



III. LINEAR PERTURBATIONS 



Aim of this section is to write down the linear perturbation equations for a combination of fluid components with 
general equations of state pi = Wi{p)pi and a scalar field with a gene ral potential V{(/)) and coupling l3i{(f>) to the 
fiuids. Some of this material has been already published (e.g. in [1^1^) but not at this level of generality. We choose 
the longitudinal gauge 

ds^ = a^[-(l + 2^')dT^ + (1 - 2<i>)dx,dx'], (20) 

where r is the conformal time. We will write the equations in the time variable a = log a, so that our metric is 
effectively 

ds^ = e2"[-(l + 2*)^ + (1 - 2<P)dx,dx']. (21) 
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It is well-known that in absence of anisotropic stress $ = ^, so we adopt this simplification from the start. We define 
the perturbation variables 

5^6p/p, ^ K5(j)/\/&, "Vi^-irrn:! ViVi = 6. (22) 

Kelt 

Repeated indexes mean summation and all spatial derivatives are with respect to the comoving coordinates Xi. We 
define also the dark energy mass 

= ^ , (23) 

with its dimensionless version 
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Tfl 

K = j^- ^^^'y'h ■ (24) 

Perturbing the Einstein equations and the conservation equations we obtain the Hnear perturbations below. We 
introduce the scale A: in real space we interpret as the operator — H~^V^; in Fourier space, A = Ti./k. In this 
way the following equations can be read equivalently in real or Fourier space. Note also that (3' — dl3 / da — (j)'dp/d(j). 

The perturbation equations for a generic equation of state p — w{p)p, which includes unified models like the 
Chaplygin gas or phase transition models are 

Generic fluid component (equation of state p — w{p)p) 

5' = -(u; + l)6' + 3(l + w)$'-2(l-3w)(/3v3'+/3V)-3w,p(5(l--2a;/3), (25) 

e' = -[(1 - Zw){\ - 2x(3) - w pA{w) + + + "^f \-H + 2/3j^^^A-V + (1 + w)\-^^, (26) 

n w + 1 w + 1 

where w,p = dw / dlog p, and 

A{w) ^3 + 2xf3{l-3w)/{l+w). (27) 

Note that the sound speed is = dp/dp — w + w^p. The equations for the scalar field coupled to several fiuids with 
equations of state pi = Wi{pi)pi and the metric equations, respectively, are 
Scalar fleld 

^" + (2 + ^)^' + (A-2 + m2)^-4$'a:-2yV/i<I' = (28) 



^A[l-3u;,-3u;,,p]r!,(5, + 2^/3,r!,$ + ^(l-3«;0-A'^^^ , (29) 



Metric 



_ -3A^[6x(^ + 2xip' - 2y^phv + E + 3(^;, + l)\^e,)] 

2(1 -3A2(a;2 +22/2)) ' ^•^^> 

$' = i[2(3a;(^-$) + A2^3(w, + 1)0,^!,]. (31) 

From these general equations one can derive the following equations for the three perfect fluid components, CDM 
(w ~ Q ~ Cs), baryons (w ~ Q, but non-negligible at small scales), radiation {w = 1/3). 
CDM 



5', = -e, + Z^' -2(3,^' -2(3',^, (32) 

n 



1 + ^ - 2l3,x ) 0, + A-2 ($ _ 2/3,^) , (33) 



Radiation 



= --0^ + 4$', (34) 
^7 = -^0, + ]x-H, + X-H, (35) 
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Baryons 

6'i, = -0b + 3$' - 2/3bV7' - 2/3^^ , (36) 
9't = - (^1 + ^ - 2/3fcxj 9b + cjX-^S + A-2($ ^ 2pb^) (37) 

(to the equation for one should add the standard term describing momentum exchange with photons due to 
Thomson scattering, seepj])- 
Scalar field 

^" +{2 + ^-)ip' + (A-2 + ni^ip - A^'x - 2yV/i<I' = 

Pc^ciSc + 2$) + f3bnbi6b + 2$) + ^{flcP', + ^bPh) , (38) 

X 

Let us derive now the Newtonian hmit (small scales, A ^ 1). The gravitational potential is 

* = -^A2(^r!,,5, + 6x(^ + 2V-2yV/iv), (39) 
$' = ixLp~<^. (40) 



Inserting ijifijl in l(38jl we obtain 

if" + (2 + ^)^' + (^(A~2 + ,^^2 _ ^2a; - ^ ^ f7,/3,') + $(4a; - 2j/V/i - 2 ^ VL^P,) = ^ , (41) 

where the sum is on the coupled components (here baryons and dark matter). It is interesting to observe that the 

terms in /?■ contribute to the equation as effective masses. If for instance Pi = /JgeV^^'"^^'^ then we can define a 
"coupling mass" 

ml = ^ = 2n,m . (42) 

X 

The Newtonian limit in Eq. lllT|l amounts to neglecting the derivatives of (p (because we can average out the rapid 
oscillations of ip) and the metric potential <& (which is proportional to A^). We can neglect also the term 12xLp since 
|a:;| < 1 is much smaller that A~^. Finally, we neglect here also and to^, with respect to A~^ since otherwise the 
dark energy would cluster on astrophysical scales and would reduce to a form of massive dark matter js^ . In section 
v., however, we remove this approximation. Prom Ij41ll finally we obtain 

p^ \\Pc^cSc + Pb^bSb). (43) 

Since cp is of order A^, Eq. Ij39|l reduces to the usual Poisson equation (hereafter we neglect radiation) 

^ = -^x''{nbSb + n,s,) . (44) 

Now, if we substitute in Ij33ll we can define a new potential acting on dark matter 

= $ - 2/3,(p = -^X^nbSbil + pbPc) - ^x^^M^ + pc) ■ (45) 
In real space, this equation becomes 

V^^c 'inGbcPbSb + AttGccPc&c , (46) 
where I restore the gravitational constant and define 

= G7y , 7y = 1 + 4A/3j/3 , (47) 
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so that Gbc — G(l + 4/3b/3c/3) and Gcc — G{1 + 4/3^/3). Analogous equations hold for the baryon force equation Ij.'-iTH . 
Therefore the Newtonian linear equations for dark matter and baryons in coupled dark energy are 























V2$b 


= 4:7rGbb. 



n 

H 



(48) 

2(3,x ] 0c - n-^V^<^c , (49) 



(50) 

2(3bx]eb-n-^w^<i>b, (51) 



inGccPcSc , (52) 
AirGbcPcSc ■ (53) 



The f3'(p terms in the S' equations have been dropped because f is of order A^. Deriving the S'^ equations we obtain 

^- + {^ + Jc~ 2'^^'') ~ ^(^cc^cf^c + IbcSb^b) = 0, (54) 



and similarly for S'^ 



^fc + ( 1 + ^ - "^/^bx] 5'b - ^{-ibc5c^c + IbbSb^b) = 0. (55) 



These equations generahze previous results because are vaHd also for non constant (3 (provided ^ •^^^)- 

It is clear that since baryons and dark matter obey different equations, they will develop a bias already at the 
linear level. A simple result can be obtained in the case in which one component dominates. Assuming Q.b ^ ric, in 
fact, the baryon solution will be forced by the dominating CDM to follow asymptotically its evolution. Putting then 
5c ~ m(a)da _ -^i^ji ^ — const. wc obtain the coupled equations 

m' + + [ 1 + ^ - 2f3cX j m- -'jcA = , 

m' + + h + — - 2PbX ] m - —jbc^c = , (56) 

from which by subtraction 



n J 2b 



+ 4(/?c - f3b)xm ■ ^ ^ 

Notice that all terms on the right hand side are in general function of time. This shows that a linear bias of gravitational 
nature develops whenever /?c 7^ Pb- This bias extends to all Newtonian scales and therefore is distinguishable from 
the hydrodynamical or non-linear bias that takes place in collapsed objects. 

Further insight can be gained when these equations have constant coefficients, i.e. when and flcb are constant 
(in [T^ we denoted these cases as stationary solutions) . As shown in 0| this is reaHzed on the critical points of 
a field governed by an exponential potential and a constant coupling f}. In this case is convenient to define a total 
equation of state 

We = Ptot/Ptot ^ - y"^, (58) 

instead of — p^/p^. Neglecting the baryons the relation is simply 

We = (1 - ^c)w,p, 

SO that now 

^ = -^[l + 3u'e]. (59) 

The scale factor in this case grows as a ^ t-p/(i-p) ^ fP where p — 2/[3(we + !)]■ The solutions are Sc — a™± and 
5b = ba™^ where b is as in (|57ll and, again neglecting the baryons, i.e. for fif, ^ f2c [3^ 

m± = i(-l + 3u;e+4/3ca:± A), (60) 
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where = 2A^cc^c + ( — 1 + 3we + ^PcxY ■ In this case then m and b are constant. The scalar field solution is 

LP « X^IScSc^c = Hlk-^a^^'P~^^/^l3c5c^c- (61) 

For small wavelengths tp (which here is proportional to ^p^j ) is always much smaller than i5c, at the present time, 
unless of course /3c is exceedingly large. It is interesting to observe that (/? could outgrow the matter perturbations in 
the future in an accelerated epoch, i.e. if p > 1 and if fie does not vanish. 

IV. SYNCHRONOUS GAUGE 

Since most Boltzmann codes in CMB are implemented via the synchronous gauge we give here the relevant equations 
in this gauge. 

Generic fluid component (equation of state p — w{p)p ) 



S' = -{w + l)e - -{w + l)h' - 2(1 - 3w)il3ip' + p'p) - 3w,p(5(l - 2x[3) , (62) 



(1 - 3iz;)(l - 2Px) - w^pA{w) + % 

n 



^^ + ^A-^ + 2/3^A-V, (63) 
w + 1 w + I 



Scalar fleld 



^" ^ ^ ^ ^ l'^'"' ^ 2^'^'-^2<^ = H - 3w« - iw^,p]nA + ^ r!,(l - 3w^)^(p , (64) 

Metric 

h' = 2\-^ri + 3^5ini + 6ip'x-6pfiy'^ip, 

V = ^X^^{w^ + l)nA + 5^x , (65) 

TT/ 

h" = -(1 + — )/i'-2[l2(y3'.T + 6Ai/iyV] -3^(l + 3iy + 3wp)5,rJi. 



Again in view of CMB applications, it is useful to detail the adiabatic initial conditions. The condition of zero 
entropy perturbations is 

5S - ^ = , (66) 

^s' = fT^V-fr^V-o. (67) 

For the scalar field 



1 + Wi J \1 + W-j 



h _ (p'<l>' + (pa'^H ^¥4, _ xLp' - ipy^pfi 

so that applying Eqs. (|6filfi7|l to the scalar field and the other components we obtain the initial conditions as 

^ xH8y^^ih6c + 2xh' + AxS',) 
^ A{x^ +y^fifiX^{6x- (3n,)) ' ^ ' 

, ^ xjAx'b, - y^iifiX'{-24x5e + 8^iflyH, + ipS^fl, + 2xh' + 4x5'^)) 

A{x^+y^fihy{6x-(3n,)) ■ ^ ^ 

In a radiation dominated era in which flc ^ and on super-horizon scales (A ^ 1), these become 

xjiy^^ifiSc + xh' + 2x6',) 

^12xS, + i^ihy^S, + xh' + 2x6', 
^ = • (72) 
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Inserting S'^ from Eq. (jf)2ll and putting initially 6'c = 0, we can further simplify 

H> = —f ^c, 73 

V. A MASSIVE DARK ENERGY FIELD. 

Here we take a digression to consider the two effective masses of the dark energy field, previously neglected. If is 
not much larger than rh^ = + , Eq. in Fourier space becomes (in this section we assume the dark energy 
is coupled to a single matter component, subscript m, or, equivalently, that has a universal coupling to all fields) 

if « Y{k)X''pnrnSm, (75) 

where 

where m = rhH. If we substitute in Ij33|l we see that the effective potential is (neglecting the baryons) 

^^^^X^aMl + ^P^Y{k)]. (77) 
Now, let us write down the density contrast for a particle of mass AIq located at the origin in empty space: 

n X - - Pm _ n^M{<t>) 
Pent sH'^a 

where 

M(0) = Moe- ^ / '^''^ , (79) 
where we used Eq. (EJ. It turns out then that the potential originated by a dark matter particle is 



which, upon inverse Fourier transform 

(2^) 



$(r) = I e'^ ^'^d^k, (81) 



becomes the Yukawa potential 



where r = ax is the physical coordinate. 

It is useful, in view of application to A^-body simulations, to write down explicitely the acceleration on particles. 
Taking Eq. Ij33|l and using the definition of peculiar velocity Vp — adx/dt in terms of the velocity used in Ij22|l 



Vp — Tiv . 



we can write the acceleration equation in ordinary time dt — adr as 



d^ 

vp^, = -{l~2l3x)Hvp,,^ —. (83) 
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If we define 



1 H --(1 + mr)e 



we obtain the force on a dark matter particle 



d$ _ GY{r)M{(j)) 
drj 



(84) 



(85) 



In Eq. H8;-ill the three effects of the couphng appear clearly: the mass M depends on the time evolution of 4>; the 
gravitational potential acquires the Yukawa correction; and the mass variation of the test particle induces an extra 
friction —2/3x. 

It is straightforward to generaHze to the acceleration of a test particle of type t due by a distribution of several 
particles of species s at distances r^t (dropping the p subscript): 



Vt=-il-2(3tx)Hwt~J2 



GYst{rst)Ms{Ps,4>) 



(86) 



where 



GYst{r)^G 



(1 + mr)e 



(87) 



In these equations have been applied to A^-body simulations in the limit of constant (3 and m — > 0. 

In practice, in any dark energy model one expects the mass scale to be much larger than the galaxy cluster scale 
K IMpc, in order to prevent clustering, so r <C 1/m for all astrophysical scales. However, a non- vanishing mass m 
can have some interesting effects. 

First, it is to be observed that the background dynamics depends on the potential and its first derivative only, while 
the mass depends on the second derivative. It is then possible to build viable dark energy models that accelerate 
the expansion but whose mass scale 1/m is between a few Megaparsecs and H^"^ = 3000 Mpc//i; in this case the 
perturbations of the scalar field would be directly observable through, e.g., weak lensing. 

Second, while the mass depends exclusively on the potential F(0), the coupling mass 



(88) 



depends on the matter content flm ■ In a inhomogeneous background, we can expect that ml will be proportional to 
Pm- This raises an interesting question, recently posed in the "chameleon" model of Ref. |23| and also in the context 
of a- varying models [2^: is it possible to have a scalar field with a large mass near a massive body like the Earth and 
a very low one in space ? This would open the possibility that scalar gravity escapes detection on Earth laboratories, 
where the Yukawa term would be exponentially suppressed if 1/m is on the submiUimetric scale, even if (3 were of 
order unity (here we neglect any possible upper bound from cosmology). In it has been hypothesized that a 
model with a constant P do in fact contain a density-dependent mass but our calculations show clearly that /?' is 
a necessary condition. This possibility will be discussed in another paper. 



VI. SECOND ORDER EQUATIONS 

Here we extend the previous calculations to second order in the Newtonian regime and in the non-relativistic limit. 
This will allow to evaluate higher order moments of the gravitational clustering. Higher order perturbation equations 
in a varying dark matter mass scenario and in scalar-tensor theories which may be reduced to particular cases of the 
present model have been studied in 

Here for simplification we consider a single matter fiuid with w = 0. In the non-relativistic limit the quadrivelocity 
remains first order 

(89) 
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The general conservation equations at second order are then (w^ — ViV^) 

6' + V,(l + S)vi = 3(1 + 5)$' - 2/3(/?'(l + S)+ 6$$' - v.V.i^' - 3$) - H^v^{l - 2l3x) 



-2/3V(l + <5+^)-^(/3"-/?'^), 
1 + ^ - 2f3x^ V, 



(l + (5-2*-2$) 



-V,($ - 2/3(^) + 2^-^ - VjVjV, + t;,(2$' + (90) 



^2 'V ^2 



where the second line of each equation contains the terms from the variation of /3(0). The scalar field equation is 

b" + (2 + ^)^'~2yV/i*](l-2*) + 



2yVV(/2 + Ai/st^) + (^'(3$' + - $'.x(3 + 6$ - G^-) - ^''a;(l - 4*) - 

^an[2^'(l - 2*) + (5 - H^v^] 



(1 + 2$) — (V.^)V,(a>-v|/) ^ (91) 



^r!„/3'(i + s,) + ^n^ifi" - 13'^) 

X IX'^ 0' 

where /a is defined in Eq. II15|I . 

Let us now derive the Newtonian limit. We can use the metric equations at first order, in particular the relation 
$ = 5" . As before, $ and <^ are of order (with respect to (5) and we neglect the time derivatives of We obtain 
from ((HDl 

J^ + V,(l + ^)w, = 2v,V,$-ij2i;2(l-2/3a;), 
l+'^-2(3x^v, = -:^yr{^-2l3^)- j^VjV.v,. (92) 

Furthermore, in the first equation we can neglect the term H^v'^{l — 2(3x) because of the non-relativistic approximation 
and the term 2viVi^ because both Vi and Vi$ are of order A at first order. Finally, in the second equation we can 
approximate 1 + <5 « 1 at the denominator in the last term, since it gives a third order correction. Therefore we are 
left with 

5'^ + V,il + 5)v, = 0, (93) 
l + ^-2/?x^«, = --Lv,($-2/3(^)-i-,V,t;, . (94) 

Applying the same approximations, for the scalar field we obtain 

^" + (2 + ^)(^' - 2t/V/i* - 2/3^!™$ 

-^^-4$'x + 2yV/2<P + 22/V'/3'/'' - pn,nS, (95) 

which neglecting the time derivatives of ip reduces to the non-linear Klein-Gordon equation 

VV - "^V - <^f^ = -l3Vlra5n'^ , (96) 
where we defined the non-linear correction 

a = 2l/^i^f3H^ . (97) 

If m and a are negligible with respect to the lenght scale A then we see that the non-linear conservation equations in 
coupled dark energy coincide with the usual non-linear Newtonian perturbation equations with an effective potential 
$ and a correction in the Euler equation due to the time variation of the dark matter mass : 

S' + Vi{l + d)v, = 0, (98) 

= ATrG„„nPmS. (100) 



1~C^ 1 

v', + {I + — - 2^x)v, + VjV = (99) 
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If atfP' is negligible but mip is not then one should use the Yukawa correction to Gmm (in this case one should also 
assume that the coefficients of the terms in <&, in H9(ill are smaller than m^). As it has been shown, Eqs. II98I1()()|I 
are valid also when /3 is a function of 0. 



VII. CONCLUSIONS 



This paper is meant to set up the formalism for future work on non-linear properties of dark energy, with an 
emphasis on its coupling to matter. Along this survey we found several results that we summarize here. 

a) We derived the general equations for the linear pertubation growth for a general dark energy potential and a 
general species-dependent interaction with matter. This generalize previous work. 

6) We derived an analytical relation between the bias induced by a species-dependent coupling and the growth 
exponent of the linear perturbations, as well as their values in term of the fundamental parameters in the case of 
exponential potential. 

c) We discussed the Yukawa correction to the gravitational potential due to dark energy interaction and we found 
that a density-dependent effective dark energy mass arises only if the coupling is non-constant. The consequence of 
this effect on equivalence principle experiments will be discussed in another paper. 

d) We derived the second-order perturbation equations in coupled dark energy and their Newtonian limit. We 
showed that the coupling introduces three corrections to the standard Newtonian fluid equations, one proportional 
to the velocity and the others which can be absorbed in the gravitational potential. These equations wil be used in a 
subsequent paper to derive the large scale skewness of coupled dark energy. 
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